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I. INTRODUCTION 


The nonlocal theory of fluids has been introduced ty 
Eringen in [1] where both localized equations of motion and a 
constitutive theory for nonlocal non heat conducting  stokesian 
fluid are formulated.  Nonlocality is introduced into the theo- 
ry in two ways. First, classical field quantities - free en- 
ergy density Ww, entropy density y], and stress tensor Tu 
- become functionals defined over the motion of the body X 
and the temperature distribution 68, and over their time and 
space derivatives (e.g. V, grad V, È etc.). They are also 
functions of the latter quantities taken at a given point 4X. 


Second, some new expressions, so called localization residuals, 
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appear in the equations of motion. The residuals are functionals 
and functions of the same character as those of the fields. They 
describe a direct long-range transfer of momentum, moment of mo- 
nentum, energy and entropy between various parts of the body. It 
is renarkable in Eringen's theory that if residuals are indepen- 
dent of € (and 0, 6, and (grad VJ ), then the body force 
residual f vanishes identically. This result, however, arouses 
two doubts. One, the assumption that residuals do not depend on 
velocity (nor on other rates) is inconsistent with the princi- 
ple of equipresence, cf. e.g. [2]. Two, the vanishing of the 
body force residual seems to be a strong restriction to the theo- 
ry since one can easily imagine a medium with non zero nonlocal 
force that could generate energy and entropy effects. 

This paper aims at presenting a nonlocal fluid where the 
long-range force does not generally vanish. It is assumed that 
the fluid is locally ideal  (incompressible, unviscous, nonpolar) 
and non heat conducting. The nonlocal effects consist in long- 
-range transfer and production of momentum, internal energy and 
entropy, and in general they depend on the distance between inter- 
acting particles, on their relative velocity and their tempera- 
tures. It is shown that even such simple assumptions lead to sig- 
nificant modifications of velocity profiles for a flow of such 
a nonlocal fluid. 


II. EQUATIONS QF MOTION 


* 


The basic equations of motion together with the Clausius- 
Duhem inequality for a nonlocal non heat conducting ideal stokes- 
ian fluid take the form (cf. [1] ): 


ev -"ęf + gradp = ef, (1) 
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ei - gh = gh - tev, 
a (1) 
4 Ch + 4:8) > gl b - i Mn i fV». 
The notation is adopted from (1) . So we have: 9 = mass dengi- 
ty, f = external body force density, p = hydrostatig pres- 
sure, & = internal energy density, h = external energy Source, 


and 7 = nonlocal body force density or momentum residual, 1 = 
= nonlocal body couple density or the residual of moment of momen- 
tum, fh = nonlocal energy production or energy residual, b = 
nonlocal entropy production or entropy residual. The body ® occu- 
pies the region © of 3-dimensional Euclidean point space Ez. 
The fields '"V(X) and 9(%), Še Y, designate the velocity and 
the fémperature of a particle occupying the position X, respect- 
ively. We assume that all the fields and residuals are sufficient- 
ly regular to possess time and space derivatives of desired order. 
A superposed dot stands for the material time derivative. Iwo 
bars above a letter signify a tensor quantity. 

The localization residuals obey the following identity, [1]: 


|: | BRT 3 . dw) = 0, (2) 
which ascertains that the global nonlocal production of momentum, 
moment of momentum, energy and entropy in the whole region VAIS 
naught. In relations (1), on the right-hand side, we have ex- 
pressions describing total effects of nonlocal transfer and pro- 
duction of the quantities mentioned above: 


¢ 
in 
B 2 2 
Ly x erp AR 
À G) 
hg = h - fev, 
^ ^ ^ 1% 1 7. 
bo = b = 5 h Bg fev 
Of course, upon (1.3) we have 1, = 0 so that 1 is fully 


a F 
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determined by A 

The body is acted upon by an external body force f, ex- 
ternal pressure p acting on the boundary ?'Y', and external rate 
of supply of heat nh. The body is forced to change its spatial 
configuration and temperature distribution. At the same time the 
processes of transfer and production of energy and entropy take 
place. They are all described by the quantities (3). In order to 
determine the response of the body to a given external action one 
must solve the equations (1) together with a set of constitu- 
tive equations which express the quantities (3) as functionals 
of the independent variables X and 8. This is the scope of 
the following section. 


III. CONSTITUTIVE EQUATIONS 


As was mentioned in the Introduction we now come to consid- 
er a medium that is locally a perfect non heat conducting fluid. 
For such a fluid we have: 


t=-pl, 8-0, + =4¥(0), = 48), (4) 


where T is the identity tensor and 4 is a heat flux vector. 
For the nonlocal residuals (3) we adopt the constitutive assump- 
tion that they are functionals of X^, V(X’) and 9(X*) over 
the wnole region ' and functions of %, V(%) and 8(%X) ata 
certain point: X. Having introduced the abbreviation 


AG) = |3%, 0%), 8(3)) (5) 
one can write the following general expressions for residuals: 


f(x) = Žane); ac] , 


i" 


Ry GO XY. LA); ^], (6) 


BD = B[AGE); AG , 
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where the bracket stands for the functional dependence on A(X’), 
x'e V, and for the function dependence on A(X). 

The constitutive equations (6) must be subject to sever- 
al conditions of general character, cf. [2] and [3). They guaran- 
tee proper behaviour of the functionals $, X, and m so as not 
to violate the fundamental principles and laws of physics. Below 
only those conditions are cited which will be referred to later 
on. The constitutive postulates state that the form of the consti- 
tutive functionals must 


(a) not violate the Clausius-Duhem inequality (1.5) (the 
entropy principle); 

(b) be invariant with respect to arbitrary rigid motions (the 
principle of objectivity); 

(c) provide the vanishing of the global effect (2); 

(d) tend to the form of a local theory in the classical contin- 
uum limit (the principle of correspondence). 


Here we add one more postulate which states that 


(e) nonlocal interactions are additive with respect to sub- 
bodies. 


The last condition is not of such general and common character. 
It follows from the assumption that the nonloca?) interactions 
are two-body interactions. One can equally well derive a nonlocal 
theory without taking that condition into account. Now let us 
consider the consequences of the above axioms. 

The condition (e) together with the natural hypothesis 
of continuity and boundedness of functionals ¥, %, and B, 
upon the Friedman-Katz theorem [4], leads to the following inte- 
gral representation of (6): 


fa) -= | serge NOD) dV(3*), 
Ry = [Fot AGO, AGO) dv(&?), (7) 


BD = | BoC ACR), ACR))- aVR’), 
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PI 


where the generating functions or influence functions x hg 
and Ja characterize the nonlocal properties of the material 
of the body while the integral represents the dependence on the 
shape and volume of the body. 

Condition (c) after substituting (7) into (2) leads 
to the following equations: 


J] fes Avv dv = 0, 
tt 


Ji Ze F(A,A)-dv? dv = 0, 
Tr | (8) 
J| € RC A) + WCA, A) av? dv = 0, 
TY 
|| CCAM) + gdy gl ALA) Dave av = 0, 
tr 
where A' = ACK?) and dV?’ = dV(X'). The above identities do 


not restrict a class of motions of the body BB , viz. the class 
of variables A. They must be satisfied for all thermomechanical 
processes, i.e. for all motions and all temperature distributions. 
What is more, from the same nonlocal material defined by the func- 
tions 7 ho and Bg one can make up bodies of arbitrary vol- 
umes and arbitrary shapes. Then, one can force them to move in 

an arbitrary way. Therefore it follows that the relations (8) 
must be satisfied identically for both all functions A(®) and 
all regions Y. If the functions A hg and By are continu- 
ous with respect to variables A’ and A , then the identities 

(8) cannot be maintained for all functions A and all regions Y 
unless the expressions under the integrals are antisymmetric with 
respect to replacing X by X” and X” by X. After some calcu- 
lations one obtains the following formulae: 


E Assum zat CA d.d 


kx? = 0, 
(9) 


icu quad did cce f 


By = BCMA) +Ż-gr- BW - dg rp 


"ep 
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where {= k’ - &, i = V(X^) - V(X) and y and Q are antisyn- 
metric under transposition A'=< A. [he first equation repre- 
sents MNewton’s third law. The second one states that the nonlo- 
cal force acts along the relative position vector 4 . The third 
and fourth ones give the general form of. generating functions hg 
and bg: respectively. As one can see, the nonlocal energetic 
effect at a point X consists of two processes: the transfer of 
energy from a particle X?’ and the production of energy due to 
the nonlocal force interaction between X? and X. Similarly, 
the nonlocal entropy effect consists of three processes: the 
transfer of entropy from a particle x’, the production of ent- 
ropy due to the energy transfer from X” to X, and the produc- 
tion of entropy due to the dissipation of energy caused by the 
nonlocal force between X' and X. 

The next condition to be satisfied is the invariance re- 
quirement (b). Under superposed rigid body motions represented 


by 


— © = Mt)sf£ + a(t), (10) 


xi 


where ü is a proper orthogonal tensor function of time and € 
is a vector function of time, it is demanded that the functions 
8, y, D and f all be unaltered apart from orientation in the 
case of vector: 


From the quantities A and A’, that is X, V, 0, P, P 
and 8”, one can compose four scalars and only one vector: 


X: antes’, eit apu? and 4, (12) 


where x? = 4:4 and 44 = XR. Therefore, the only possible 
representation of f, % and- „By is: 


alre OG rakes 012-4. 
(13) 
Pede, W gau, 


am 
" 


482 X. J. Szymona 


B = ©(%,%, 9, 8' D = (115) 


and @, y and f are arbitrary scalar-valued functions with 
respect to (10). Of course g is symmetric and y and f are 
antisymmetric upon transformation x — -X and 8’ «— 8. One 
can easily see that the functions R and by possess appropri- 
ate symmetries and invariance properties as well. 

Let us now proceed to determine the consequence of the ther- 
modynamic admissibility. Condition (a) ascertains that every 
solution of the equations of motion (1.1 - 1.4) must identically 
satisfy the Clausius-Ouhem inequality (1.5). After substituting 
(13) to (1.5) we obtain: 


- t 9 ZĘ + m )8 + 


(14) 
+ e JC- B -a-p ip RAR) oV > 0. 
r 
This inequality must hold for arbitrary fields 8, 8 and 1. 
Consider first a rigid motion of the body ® . For such a motion 


we have X = 0, and so the expression under the integral vanishes. 
Thus, we get the classical relation between entropy and free en- 


ergy: 


4 = - SĘ. (15) 


Now the inequality (14) is reduced to the expression with the 
integral only. Upon using the same arguments as those exploited 
in the consideration of condition (c) we conclude that the 
inequality (14) must be satisfied locally, that is the expres- 
sion under the integral must be non-negative. Having taken into 
account the Symmetry properties of Q, y and Q we obtain 
finally 


-a-p + EGre Pre > ode > 0. (16) 


which is the condition for thermodynamical admissibility of gen- 
erating functions (so % and (5. In particular, if one as- 
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sumes that within the body there is a nonlocal force only (9 # 0), 
other interactions being naught (y= (= 0), one gets 


Ta ef BB) at 20. (17) 


This inequality states that the nonlocal force is dissipative so 
one may call it the nonlocal viscosity. Of course, the nonlocal 
viscosity force gives rise both to nonlocal production of inter- 
nal energy and to nonlocal production of entropy: 


Do. rd: et oa 
(18) 
Doa lg. XX et 


Both of them are positive and each of them is the same at both in- 
teracting points X and X'. 

Finally, we must analyse the principle of correspondence 
(d). Nonlocality can be roughly characterized by a parameter a 
which stands for a range (or a mean range or an effective range) 
of the nonlocal interaction. Similarly, every problem of motion 
is characterized by a parameter ł being a typical length in 
this problem. The specific nonlocal effects can only be detected 
in the range of ł/a ~1. When ł/a »1, then the nonlocal theo- 
ry should not differ from its local counterpart. In the limit 
a/ł —= 0 the equation (1.2) should turn into the well known 
Navier-Stokes equation. Thus, one must have the following asymp- 
totic behaviour of the generating function Pi: 


e? [oca e OEM) Spe mca) Ave, 0» 
YP 
where m(@) is a (local) viscosity coefficient. 

Similar conditions must be imposed on the generating func- 
tions of energy and entropy residuals. In a classical theory limit 
one should obtain a classical expression for the energy production 
due to a heat flow, e.g. -xæ div B, where » is a heat conduc- 
tivity coefficient and § = grad 8. For the entropy flux the 
situation is similar so we will not pay more attention to it. 

It is obvious that the ranges of y and @ may in general be 
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different from that of Q. 
Summing up the above considerations we conclude that with- 


in a nonlocal medium where the generating functions of residuals 
depend on the distance between particles, on their relative ve- 
locity and their tenperatures there may occur a priori three 
groups of effects: 


(1) Any two particles interact with each other with a long- 
-range force. This interaction gives rise to dissipation 
of energy and to production of entropy. 

GE Any two particles exchange energy in a long-range way. 
The energy flows from the particle of higher temperature 
to the particle of lower temperature. That flow is accom- 
panied by a production of entropy. 

(III) Any two particles exchange entropy by a long-range ent- 
ropy transfer. The entropy can flow both from the parti- 
cle of higher temperature to the particle of lower tem- 
perature and the other way round. This flow of entropy 
can only appear when the productions of entropy named 
above under (I) and (II) take place and its absolute 
value cannot exceed them. 


The processes presented above are described by the generating 
functions Qf, 4 and Q which satisfy condition (15) and 
appropriate conditions resulting from postulate (d). In the next 
section we exploit them to solve an example of a flow of nonlocal 
fluid. 


IV. POISEUILLE FLOW THROUGH A CAPILLARY 


Here the results of calculations af velocity profiles for 
a stationary laminar flow of a nonlocal fluid through a circular 
capillary are presented. In order to perform the calculations one 
must choose a shape of generating function qQX. Let it be, with- 
in the first approximation limit, linear in X and independent 


of temperature: 


On a Perfect Fluid with Noniocal Interactions 485 


F(R) = e? | a IBGO = PD) away, (20) 
Y 


Conditions (16) and (19) lead to the following formulae: 


NOTY PR 05 
(21) 


6 
[iate = gm. 
0 - 


The flow is induced by a constant pressure gradient and an inter- 
action with the boundary. Remaining within the frames of a non- 
local velocity-dependent model we assume that the interaction 
with the boundary, viz. the external region cv ; is also nonlo- 
cal and has the same form as the interaction within the fluid. 
Thus, we choose 


foy sgh seyn M Un) OP qva, — (22) 
cr > 
where V(X) is a definite velocity field outside the capillary. 
Further we put V = 0. Function 8 satisfies condition (21.1) 
and condition (21.2) with some different nonlocal range b and 
viscosity coefficient Y . Both functions A and B satisfy 
additionally the postulate of neighbourhood, cf. [5], according 
to which the influence on a particle X caused by other parti- 
cles decreases with the distance between particles. Thus, the 
functions A and B have maximum for y= 0 and they tend to 
zero quickly enough with 4 — sa. As it is suggested by the 
results of other calculations, cf. [5], the specific form of the 
generating functions is not as crucial as the above general prop- 
erties. We shall not go into details of calculations as this is 
not the aim of this paper. After substituting (20) and (22) 
into (1.2) and performing transformations and integrations one 
obtains only one equation 
4 
USEMANAĆ Ferre SPURS rely 404! ):ów'. (23) 
0 


Here § and % are some definite functions connected with A 
and B. The kernel % is posítive-definite and nonsymmetric. 
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wim) = ¢ M /(9g n2)) v(r) is the dinensionless velocity, w = 


r/h 


is the dimensionless radial coordinate in the appropriate 


cylindrical systen of coordinates, - 99 = ?p/9z is the pressure 
gradient in the direction of the flow, h is the capillary radius, 


a 
* 


The results of the calculations 
in the form of velocity profiles 
are shown in the figures as func- 
tions of three parameters: the 
capillary radius h/a, the range 
of interaction with the boundary 
b/a, and the intensity of the 
interaction with the boundary v/u. 
As can be seen in the figures, 
for wide capillaries (h/a> 1) 

the profile lies near the para- 
bolic curve of the classical 
flow. The interaction with the 
boundary causes a strong resist- 
ance to the flow near the wall 
within a layer of thickness of 
about the range of interaction 
with the boundary. The effect is 
very conspicuous for short and 
strong interaction with the wall 
(b/a = .5, V/M = 4). Of course 
the deviations of the profile 
give rise to changes of the 
volume rate of flow as compared 
with that of the local fluid. 

The estimation of the form of 
influence functions A and B 
must be obtained upon considera- 
tion of more complicated prob- 
lems (e.g. time-depemdent flows, 
turbulent flows, flows through 
capillaries of a finite length 

or of non circular section, etc.) 
and upon comparison with experi- 
ment. 
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STRESZCZENIE 


Przedstawiono model cieczy nielokalnej. W cieczy występują 
zależące od prędkości i temperatury zjawiska siły dalekiego za- 
si$gu, dalekozasiégowego przekazu energii i dalekozasięgowego 
przekazu entropii. Otrzymano warunki dopuszczalności dla funkcji 
tworzących residuów. Profile prędkości dla przepływu Poiseuilla 
przez kapilarę wykazują różnice w stosunku do profilów klasycz- 
nych. 
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IpuBOZKTCA MOAGJP HENOKAJIBHOŃ XMAKOCTM, B XMAKOCTM NMpoABIA- 
OTCA, S38BMCHUMe OT CKODOCTM M TeMHepaTypH, 3QQEKTH CHA AanbHero 
AelioTBMR, AANBHCHACTAHNKOHHNE NepeĄAUK 3Heprus M 3ATPONKK, Nony- 
JgeHH JCIOBKA NPKMCHMMOCTH ANA QYHKUKK OÓpa3yDMKX pe3UAYOB. Ipo- 
dunn ckopocte ana noToka IlyagefilA uepe3 KaNKIIApP NpOABIACT pas- 
NMYKE NO OTHOUEHMD K KJIACCKHUECKUM HPOÓKIAM, 
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